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Abstract 

Nuclear quadrupole resonance measurements have shown evidences that the heavy fermion 
compound CeRhlns exhibits a coexistent phase with commensurate antiferromagnetism and li-wave 
superconductivity. In order to clarify the nature of the spin-excitations in the coexistent phase, 
we have applied the RPA method to an itinerant model, where the effective interaction is given 
by two mean-field terms of commensurate antiferromagnetism and d-wave superconductivity. It 
is shown that, around the transition line between the antiferromagnetic and the coexistent states, 
a low-energy incommensurate spin-excitation is found to develop due to Fermi surface nesting. 
This feature reminds of the switching of magnetic ordering wave vector observed in the neutron 
diffraction. Further, we also calculate spin relaxation rate, which gives a reasonable explanation 
of the temperature dependence of NQR relaxation rate in the system with the coexistent ground 
state. 

Keywords: heavy fermion superconductor, antiferromagnetism, d-wave superconductivity, coexistent phase, 
spin-resonance, random phase approximation . . . 
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I. INTRODUCTION 



The strongly correlated electron system shows a rich phase diagram around an antifer- 
romagnetic quantum critical point. The phase diagram consists of a paramagnetic phase 
showing a non-Fermi liquid behavior, an antiferromagnetic phase, a d-wave superconduct- 
ing phase, and a coexistent phase of antiferromagnetism and superconductivity. Actual 
compounds showing the coexistent phase are UPd2Al;^, Celn^, CePd2Si2^, CeRh2Si2^'^, 
CeCu2Si2^i^, CeRhlng^i^, and BaFe2(Asi_a;Pa;)sr^- Among these compounds, UPd2Al3 with 
more than two f-electrons is believed that two localized f-electrons contribute to the an- 
tiferromagnetism, while remaining f-electrons with a large mass form Cooper pair in the 
coexistent phase. On the other hand, for other compounds, the same f-electron plays both 
roles of antiferromagnetism and superconductivity in correlated metal. 

Since the discovery of unconventional superconductivity of CeRhlus in 2000, this com- 
pound has been investigated intensively by many measurements under pressure. In the 
ambient pressure, the compound shows the antiferromagnetic phase around 3.8 K, in which 
the staggered moments of 0.6 yU^/Ce parallel in the a-b plane align with an ordering wave 
vector (vr, vr, 0.67r)ii""— . The antiferromagnetic transition temperature decreases with apply- 
ing pressure. Above 2 GPa, unconventional superconductivity has been observed below 2.1 
K. 

In the superconducting state of other compound with the same crystal structure 
CeCoIns, the so-called resonance peak has been observed by the inelastic neutron scattering 
experiments^. A coexistent phase of antiferromagnetism and unconventional superconduc- 
tivity is also suggested by the nuclear-quadrupole-resonance (NQR) measurement^^. 

Recently, the magnetic neutron diffraction^^ and NQR measurements^ have been carried 
out under the pressure to study the compound around the coexistent phase. The magnetic 
neutron diffraction measurement has showed a change of magnetic ordering wave vector from 
the original one at the ambient pressure to (vr, vr, O.Svr) in the low temperature region above 
1.5 GPa. This implies that an additional spin mode different from the original one at the 
ambient pressure exists in low energy spectra. Actually, the evidence of this from different 
measurement is an upturn of NQR relaxation rate devided by temperature 1 / (TiT) in the 
low temperature limit at 1.3 GPa, where the compound will be in the original magnetic 
state. It should be noted that the NQR spectrum in the coexistent phase is reproduced well 
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by an assumption of the commensurate magnetic ordering. Further, the value of 1/ (TiT) 
in the low temperature limit in the coexistent phase at 1.5 GPa is two orders larger than 
that in the superconducting phase at 2.3 GPa, where the quantity becomes temperature 
independent in the low temperature limit of both phases. 

The main motivation of our study is to clarify the low energy spin excitation around 
the coexistent phase of antiferromagnetism and unconventional superconductivity. In the 
antiferromagnetic phase, the antiferromagnetic spin wave appears as a Goldstone mode in 
the rotationally symmetric spin system. On the other hand, the resonance mode is observed 
around an antiferromagnetic ordering wave vector in unconventional superconductors. In 
the coexistent phase, both collective modes are expected. Further, recalling the analysis 
of NQR spectrum in the coexistent phase, the reason of the stabilization of the coexistent 
phase with the commensurate magnetic ordering is required. 

This paper is organized as follows. In section |lTl we introduce the mean- field Hamiltonian 
with both of the AF exchange and the BCS interaction terms and present the self-consistent 
solution in the coexisting state of the AF and the SC order. In section IIIIl we carry out 
the calculations of the dynamical spin susceptibility in the presence of the AF and the SC 
long-range order. Section HVl shows the NMR relaxation rate 1/Ti in the coexistence(CO) 
phase. Section |V] is a conclusion. We put some technical details in appendixes. 

II. SELF-CONSITENT MEAN-FIELD EQUATIONS IN THE COEXISTENCE OF 
ANTIFERROMAGNETISM AND SUPERCONDUCTIVITY 

In this section, we introduce an effective Hamiltonian describing the CO phase of an- 
tiferromagnetism and d-wave superconductivity in a quasi-two-dimentional system, which 
will be a simplified version of CeRhlus, and show the one electronic state in the CO phase. 
Recalling that the staggered moment is observed in the antiferromagnetic state of the com- 
pound, the antiferromagnetic phase is a metallic state with the Fermi surface. Further, it 
has been shown by a detail analysis of the inelastic neutron scattering data that other com- 
pound CeCoIns with the same crystal structure shows dx'i-y^- superconductivity at 2.3 K 
in the ambient pressure^. Since the superconducting phase of CeRhlns appears around the 
antiferromagnetic phase and the transition temperature is quite close to that of CeCoIns, the 
(i2;2_j;2 -symmetry of superconductivity in CeRhlus is plausible. Assuming the single orbital 
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for electronic state in the unit cell of the paramagnetic phase for simplicity, the effective 
Hamiltonian of the system is given by 

H = Hf, + Hsc + Hjng (1) 

where H^, H^c, and H^g are the kinetic energy, the interaction of rf-wave pairing, and the 
antiferromagnetic Heisenberg-type spin exchange terms, respectively. 

Applying the tight-binding approximation, the non-interacting kinetic term is given 

as 

Hk = ^^Skcl^Cka, (2) 
k cr 

with the dispersion relation measured from the chemical potential fi 

£k = — 2t(cos kx + cos ky) — At' cos k^ cos ky — /x, (3) 

where Ck,(7 is an annihilation operator of a quasi-particle with momentum k and spin a. 

In order to describe superconductivity, the BCS Hamiltonian is used within the mean-field 
theory 

Hsc = -^5^5^{^k,k'(cLcLk_.)c_k'-aCkv (4) 

k,k' cr 

+ Vk,k'(C-k'-<xCk'<x)cLclk_^}, (5) 

where the pairing interaction Vk^k' is assumed to take following form 

Vk,k' = V(l)k(l)k' (6) 

with a basis function of dx2_y2 wave gap function 0k = cosk^ — cos ky. In the following, 
(■ ■ ■ ) is the thermal average with the effective Hamiltonian. 

For the antiferromagnetism, the mean-field approximation is applied for the Heisenberg 
spin exchange, 

which is characterized by the staggered magnetic moment, 

'^Q = ^EE^kX.'Ck+Q.'. (8) 

k o-cr' 



The order parameters for the superconductivity 

1 
2 



Ak = 5^ ^k,k'a(cLk_.4) = A0k (9) 



k'o- 



and the antiferromagnetism 



(10) 



are self-consistently determined by solving the equations of motion for Green's functions, 



|:G(k,r) = -<5(r)-(£k-/x)G(k,r) 



d_ 



+ ?7m,F-~Q(k,r) + AkF^t(k,^), 
F^Q(k,r) = -(ek^Q-;.)F^Q(k,r) 

+ [/m,G(k, r) + Ak+QF*Qt(k, r), 



ar 

+ ?7m,F*Qt(_k,r) + AS;G(k,r), 
|:F*Qt(k,r) = (ek+Q-/i)F*Qt(k,r) 

+ ?7m,F^t(k,^) + A*+QF^Q(k,r) 

where the Green's functions are defined as 

G(k,r) = -l^(r.[cu.(r)cL(0)]), 

(J 

^^''(k,r) = -^$^^(T.[ck+Q.(r)4(0)]), 
F^(k,r) = -^5^a(r.[ck.(r)c.k-.(0)]), 



(11) 
(12) 

(13) 

(14) 



(15) 



with Cko-(T) = e^'^Cko-e"'^^. G'(k, r) is the normal Green's function, while F^^(k, r) and 
F*(k, r) are the anomalous Green's functions associated to the antiferromagnetism and the 
superconductivity, respectively. The last Green's function F*^^(k, r) proportional to the two 
different order parameters takes a non-zero value only in the CO phase, and is the so-called 
TT-triplet pairing. 

Following the symmetry arguments given in Table [U the vr-triplet pairing is clarified 
as an antisymmetric pairing in two sublattices. Since the pair of two electrons is a wave 
function of two electrons, it is antisymmetric for the exchange of the two electrons. In 
the paramagnetic state, this antisymmetry is satisfied in the spin-singlet pairing with even 
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triplet 
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singlet 
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antisymmetric 
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triplet 


even 


antisymmetric 



TABLE I: Pairings in the coexistent phase of antiferromagnetism and superconductivity. 



parity or the spin-triplet pairing with odd parity. On the other hand, the sublattice index is 
additionally necessary in the phase with the antiferromagnetic order parameter. Therefore, 
the spin-triplet and sublattice-antisymmetric pairing with even parity and the spin-singlet 
and sublattice-antisymmetric pairing with odd parity are also possible in the coexistent 
phase. Since the parity is preserved, the spin-triplet and sublattice-antisymmetric pairing 
with even parity is induced in the present case. 

The self-consistent mean-field equations for the electron density rze , the staggered mag- 
netic moment irts and the superconducting gap function are obtained as follows.—"— 



T 



No 



(16) 



■ tanh ■ 



2Eu 



No 



2T 2Eu 



tanh ■ 



2T 



(17) 



Eu 



No^ y/t:i^ + {Um,y2E^^ 
k' 



tanh 



tanh 



E, 



k+ 



2T 
2T 



E^k,k' 



Ak^ 
2E\ciA 



(18) 



tanh ■ 



E-ki^ 



2T 2Eu'- 



Ak' , , E]^/^ 
tanh ■ 



2T 



where Yl'k is the momentum summation within the magnetic Brillouin zone. Here, the 
dispersion relations of quasi-particle in the coexistent phase is obtained by solving the secular 
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Fig. 1 (Color online): Phase diagram obtained by solving the self-consistent equations, Eq. (jl7p 
and Eq. (jlSp with a set of parameters U = 2.2, V = 0.8 and t' = —0.2 in units of the nearest 
neighber hopping t. The trajectory with open symbol shows the temperature evolution of the 
system for the fixed chemical potential /i = —0.69 and = —0.85. 

equation of four Green's functions as follows. 



are the excitation-energy in the antiferromagnetic state with = —2t{coskx + cos ky) 
and = — 4:t' cos cos ky. The detailed calculations for Eq. f lTB]) ^f l2U]) are presented in 
Appendix lAl The phase diagram of the system with a parameter set t' /t = —0.2, U/t = 2.2, 
and V/t = 0.8 is shown in Fig. 1 by solving the self-consistent mean-field equations. Each 
phase boundary in Fig. 1 is obtained with tunning T or Ue around the critical values, around 
which Eq. ( IT7|) and Eq. ( IT8|) can be linearized with respect to the small order parameters 
nis and Ak. For examples, the phase boundary between the CO and the pure SC phase can 




(19) 



where 




(20) 
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Fig. 2 (Color online): The AF and the SC order parameters and A for the chemical potential 
H = —0.69 and fj, = —0.85 in units of the nearest neighbor hopping t. At /i = —0.69, the average 
magnetization is slightly suppressed when the SC long-range order coexists with the AF. 



be tunned by the linearized equation 



U 
2iVo 



U 1 

2iVo 



: tanh 



/2 
k+ 



2T 



: tanh 



|AkP 



2T 



(21) 



where i^ki = 'f2k — At =t l^ikl- On the other hand, the phase boundary between the CO and 
the pure AF phase can be tunned by the linearized equation. 



2K 



tanh 



k+ 



2T 2Ek- 



tanh 



2T 



(22) 



where Eit_± is the excitation energy in the AF phase defined in Eq. (1201) . 

In the coexistent phase, one order parameter is affected by another order parameter. The 
effect of A ^ on the magnetic order parameter rris in the CO phase is shown in Fig. 2(a) 
with the fixed chemical potential at fi/t = —0.69, whose trajectory is also shown in the 
phase diagram. The existence of A 7^ reduces below T < Tc as shown in Fig. 2(a). On 
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Fig. 3 (Color online): (a) A dispersion relation E]^^ = t2k ~ P' ~ y ^ik {Uvrts)'^ crossing the 
Fermi level with nis = 0.132, T = 0.0347 and ^ = —0.69 in the AF state, (b) Density of states 
p{E) of the £^k_-band in the AF phase. 

the other hand, a similar figure of A in the pure SC state with /i/t = —0.85 of the phase 
diagram shows usual behavior of mean-field order parameter. Although the suppression of 
mg in the CO phase is less than 10 percent of the full magnetism in the pure AF phase, the 
SC gap opening at the Fermi surface makes a significant change into the density of states 
near the Fermi level as shown in Fig. 3(b) and Fig. 4(b). 

Fig. 3(a) and Fig. 4(a) show the dispersion of the excitation energy in the AF and the CO 
phases in the 2D Brillouin zone (BZ), respectively. For convenience, we plot only i?k- and 
— E'k- bands located near the Fermi level. In the AF phase(Fig. 3(a)), there are four Fermi 
surfaces with E\^_ = around the magnetic BZ boundary, \kr^\ + \ky\ = ir, each of which 
encompasses a hole-pocket with i?k- > 0. In the CO phase(Fig. 4(a)), the Fermi surface 
is split by the SC gap Ak, whose maximum A^""^ = 0.024 is less than one percent of the 
band width W ~ 2.5. We see that the SC gap opening causes sharp coherence peaks near 
the Fermi level, with which the spin-excitation spectrum exhibits considerable enhancement 
by a nesting effect, as shown in Fig. 5(b). We discuss the issue in the following section in 
terms of the dynamical spin susceptibility calculated with the RPA. 
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Fig. 4 (Color online): (a) An energy branch — -Ek- = ~y^k- ~^ l^kp below the Fermi level in 
the CO phase with = 0.132, A = 0.012, T = 0.0196 and = -0.69. (b) Density of states p{E) 
of the —Ei^^-hand in the CO phase. 

III. DYNAMICAL SPIN SUSCEPTIBILITY IN THE PRESENCE OF ANTIFER- 
ROMAGNETIC AND SUPERCONDUCTING LONG-RANGE ORDER 

In order to study the magnetic correlation around the coexistent phase, we calculate 
dynamical spin susceptibility. In the AF and the CO state with rris 7^ 0, there are two sites 
of different sub-lattices in the magnetic unit cell. Then, the area of magnetic Brillouin zone 
is a half of the original one, and the ordering wave vector Q becomes a reciprocal lattice 
vector. Therefore, momenta of initial and final states of relevant spin fluctuation can differ 
by the reciprocal lattice vector as the Umklapp process. In this case, the spin fluctuations 
have off-diagonal terms in momentum space representation. Since this fact is independent 
of components of spin operators, we introduce a generalized form of the fluctuation with 
a, /3 = 0, X, y, z, components'^, as follows, 

ri/T 

X"''(q,q',ifi„)= / (ire'^-M^r 5S^{t)5S^_^,{0) ) (23) 

Jo L -I 

where S^ = ^ J2k J2a '^ko-'^k+qo- is the charge operator with a momentum q. 

In a phase with a staggered moment 7^ parallel to ^-direction, the transversal com- 
ponents of q',if^n) with a, P = x,y are decomposed from the charge and longitudinal 
components. The transversal components have the contribution of spin-waves dominating 
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(a)T=0.0347t 



(b)T=0.0196t 




(c)T=0.0178t 





(d)T=0.0122t 




Fig. 5 (Color online): The momentum-dependence of the spin-excitation spectrum 
Imx^p^(q, in the limit of a; — )• after the analytic continuation iil„ ui + 16 at four different 
temperatures T =(a)0.0347, (b)0.0196, (c)0.0178 and (d)0.0122 in units of the nearest neighbor 
hopping coefficient t. The chemical potential is fixed at /i = —0.69, with which superconductivity 
develops at T = = 0.0196. The singularity at q=Q appearing at all temperature is the Goldstone 
mode due to the AF long-range order, (b) At the critical temperature T = Tc, the spin-excitation 
spectrum exhibits another incommensurate peaks at q=Q-|-(5q with (Jq/vr = (±0.08, ±0.06) and 
(±0.06, ±0.08). The momentum of the incommensurate peaks agrees with the nesting vector, as 
shown in Fig. 6. 

the low-energy spin excitation, while the longitudinal component of ^^^(q, q'^i^n) reduces 
to be massive due to the staggered moment. Especially, when the system is rotationally 
invariant, as in the present case, the spin wave becomes a Goldstone mode at the order- 
ing wave vector Q. Therefore the spin-excitation spectrum in the presence of the magnetic 
long-range order, which can be probed by the inelastic neutron scattering (INS) and the nu- 
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Fig. 6 (Color online): Nesting vector Q + 6q 

clear quadrupole resonance (NQR) experiments, is dominated by the transverse component 
^-'-(q, if2n), especially, around the superconducting instability in the AF state. In the 
RPA calculations of dynamical susceptibility in the CO phase, the transversal component of 
the RPA susceptibility XRPAi^i i^n) is written as a 4 x 4 matrix consisting of the transverse 
component of the irreducible susceptibility x"'^ {c[, with a,(3 = x,y as described in 

Appendix [B] in detail. Even within RPA, the Goldstone mode appears in both AF and CO 
phases. The equality of the divergence of transverse susceptibility with the self-consistent 
equation of the staggered moment in Eq. ( 122|) is analytically shown in the coexistent state 
of our model, as in the antiferromagnetic case^^. 

Figure 5 shows the momentum- dependence of the spin-excitation spectrum 
Imx;^p^(q, in the limit of a; —t- at four different temperatures. The chemical po- 

tential is fixed at /i = —0.69, with which a system shows transitions to the AF and the CO 
states at T/v = 0.132 and Tc = 0.0196, respectively. In Fig. 5, every temperature is chosen 
below T/v so that one can observe the Goldstone mode at Q=(7r, vr) 

lini lini Imx^PA(q, = oo. (24) 

At T = Tc = 0.0196, where d-wave superconductivity starts to coexist with the AF, we find 
sharp enhancement of spin excitation spectrum at q = Q + (5q with 5q/7r = (±0.08, ±0.06) 
and (±0.06, ±0.08). The incommensurate spin excitation at q=Q+5q turns out to be a 
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nesting contribution as shown in Fig. 6^. Such peaks appear only in the vicinity of T = Tc 
(Fig. 5(b)) and are rapidly suppressed with lowering T below the transition temperature 
(Fig. 5(c), 5(d)), because the electronic states around the Fermi level is gapped in the SC 
state. Decreasing temperature below T = the frequency of incommensurate spin reso- 
nance shifts from around zero to a finite value smaller than 2A, as shown in Fig. 7(c) and 
7(d). Recalling the spin resonance stabilizes the d-wme superconducting phases of high-Tg 
cuprates and CeCoIus, the coexistent state of our model will be stabilized by the incommen- 
surate spin resonance, which is characteristic of the coexistent phase. In the following, 
we compare this result with experimental data of CeRhlns to give a comment. The neutron 
diffraction experiment shows that the ordering wave vector approaches from an incommen- 
surate one at the ambient pressure (vr, vr, O.Gvr) to another incommensurate one (tt, tt, O.Svr) 
closer to the commensurate one at 1.48 GPar^. Similarly, it has been reported that the NQR 
spectrum also change a broad shape interpreted by assuming an incommensurate order at 
the ambient pressure to a sharp peak form explained by the commensurate order at 1.86 
GPa, where this compound is in the coexistent phase in the low temperature region^. On 
the other hand, it has been shown by the dHvA experiment that the /-electronic character 
also changes from the localized one at ambient pressure to an itinerant one above 2.3GPa 
due the Kondo effect^^. Considering these, switching of the ordering wave vector from the 
incommensurate to the commensurate one seems to be attributed to the Kondo coupling 
between the /- and conduction-electrons, which changes the shape of the Fermi surface of 
the quasi-particles. Now, the sharp incommensurate spin excitation intensity at T = Tc 
shown in Fig. 5(b) means an instability of the commensurate antiferromagnetism, though 
the order parameter is of the commensurate antiferromagnetism. Around the instability, 
the free energy of corresponding incommensurate antiferromagnetic state may be lower than 
that of the commensurate one so that an first-order phase transition from the commensurate 
to the incommensurate ones is expected in our model. Then, the incommensurate antifer- 
romagnetic state will have some similarity with the incommensurate ordered state observed 
by the neutron diffraction measurement under the pressure, because heavy quasi-particles 
around the Fermi level will be responsible for the incommensurate antiferromagnetism, where 
the itinerant picture using quasi-particles will be applicable. In addition, it is desirable to 
observe the spin resonance mode characteristic of the coexistent phaes as shown in Fig. 7. 
Finally, we give another consideration on how the coexistent phase with the commen- 



13 




Fig. 7 (Color online): The spin-excitation spectrum Imx^p^(q, w) in the q-uj plane along the 
diagonal direction = qy at four different temperatures T =(a)0.0347, (b)0.0196, (c)0.0122 and 
(d)0. 00686 in units of the nearest neighbor hopping coefficient t. The chemical potential is fixed 
at /i = —0.69, with which superconductivity develops at T = Tc = 0.0196. The strong signal at 
q=Q appearing at all temperature is the Goldstone mode due to the AF long-range order. The 
incommensurate zero-energy mode appearing at T = Tc in Fig. 5(b) and Fig. 7(b) converts to a 
strong resonance peak as temperature decreases below T = Tc (Fig. 7(c) and 7(d)). 




FIG. 8: Feynman diagram of AF*Qt(A:). The dashed line means the on-site Coulomb repulsion U. 
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surate magnetic order is stabilized. As the stability of coexistent phase has been already 
discussed, the sign of a coupling term, like 7I Akpm^ + 7AkAk+Q?Ti^ between two types of 
order parameters in the free energy, determines whether the coexistent state is favorable (- 
sign) or not (+sign)22.. Similarly, the magnitude of the vr-triplet pairing F*'°^'^(k, r) will give 
an important information, because the 7r-triplet pairing is proportional to the product of 
two different order parameters within the mean-field theory. Beyond the mean-field theory, 
a diagram shown in [8] contributes to the vr-triplet pairing as follows, 

AF'^\k) = -G(°)(k+Q,-iu;„)G(k,iw„) x 

^'^'W - + ^''(^ - A:'))F*Qt(A:'). (25) 

Since the spin fluctuation exchange term includes the spin-wave mode, the correction is 
expected to enhance the mean-field value of F^^'^{k). Through the correlation processes as 
given here, the coexistent state with the commensurate magnetic order will be stabilized. 
Actual calculation is remained future problem. 



IV. SPIN-LATTICE RELAXATION RATE IN THE COEXISTING PHASE 

The spin-lattice relaxation rate 1/Ti can be written as 

1 ^ ^ 

T^ = l'A' {S^{t)S^me'-''dt (26) 
J- 1 Jo 

where 7 is the gyromagnetic ratio of the nuclear spin I and A is the hyperfine coupling 
constant between the nuclear and the quasi-particle spins S, given by T-ihyp = ■ S. The 
nuclear Larmor frequency can be set to be zero in the limit ujqT ^ 1 where r is the 
fluctuation time scale of the hyperfine field. 

The expression in Eq. (l26l) can be rewritten in terms of the dynamical susceptibility 
through the fluctuation-dissipation theorem as — 

1 = ^ hm ^ y ^2(q)^^WSli^±!2:) (27) 

Figure 9 shows the spin-lattice relaxation rate 1/Ti calculated in systems of SC and CO 
phases with fixed chemical potentials iijt = —0.85 and —0.69, respectively. For each case, 
the trajectory in the T-rig phase diagram are shown in Fig. 1. For [ijt = —0.85, a system 
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Fig. 9 (Color online): The NMR spin-lattice relaxation rate 1/Ti as a function of temperature 
for the fixed chemical potential /x = —0.69 and —0.85 in units of the nearest neighbor hopping 
coefficient t. Figure 1 shows the trajectories in the phase diagram for the two cases, which end up 
in the coexistence (CO) and the superconducting (SC) phases at T = 0, respectively. 

enters the SC phase around Tc/t = 0.09. Therefore, the T^-behavior is due to the dispersion 
relation of quasi-particle 1/5^ + superconducting state with d^2_y2 symmetry. On 

the other hand, for fi/t = —0.69, it undergoes successive transitions at T^/t = 0.132 and 
Tc/t = 0.0196 to arrive at the CO phase. Similarly, the dispersion relation of a quasi- 
particle branch -sj + is responsible for the T^-temperature dependence below Tc for 
fi/t = —0.69, while the NQR relaxation rate shows the metalic behavior due to the existence 
of the Fermi surface in the temperature range Tc < T < T/y. For both cases, 1/Ti behaves 
as T-linear in the lowest temperature region due to the finite imaginary part 6/t = 0.001 in 
the process of analytical continuation ifi„ — )■ a; + 15. This point is discussed in the following. 
These behaviors shown in Fig. 9 are entirely consistent with NMR experimental data in the 
superconducting and coexistent phases of CeRhlus. 

The crossover temperature to 1/Ti oc T will be determined by the effective coupling 
constant of impurity scattering. In our calculation, the T-linear behavior is originated from 
a small but finite imaginary part 6/t = 0.001 of frequency in the analytic continuation 
ifln — )■ w + i^^S. On the other hand, considering the impurity scattering of quasi-particle, 
the self-energy correction will be given by Sj^p = —i6imp = — i7™jmp|%{:k'P'^(-^k- — Ew~) in 
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the present case, where Uimp and u^-' are the impurity density and the Fourier transformed 
effective impurity potential, respectively. If the impurity is included in the sample, the 
residual density of states around the nodal points of the superconducting gap is induced 
by the impurity scattering. Then, the Korringa behavior is expected from the residual 
density of states^"—. Here, considering the vertex correction by the spin fluctuation for 
the impurity scattering, the magnitude of effective impurity potential is enhanced^^. The 
crossover temperature measured in the NQR experiment is almost independent of pressure 
varing from 1.82 to 2.35GPa, where the ground state changes from the coexistent state 
to the superconducting state by increasing pressure. Thus, the contribution of the vertex 
correction by spin fluctuations will be small, because the pressure of 2.05GPa is really close 
to the critical pressure of the antiferromagnetic quantum critical point. 

V. CONCLUSION 

We have formulated the spin-excitation spectrum in the coexistent phase of antiferro- 
magnetism and (i-wave superconductivity within RPA to calculate the low energy spin exci- 
tation. It has been shown that additional low energy spin excitations with incommensurate 
momenta develop around the phase boundary between the antiferromagnetic phase and the 
coexistent phase. This implies that the corresponding incommensurate antiferromagnetic 
phases can appear between the antiferromagnetic and the coexistent phases^^i^. We have 
suggested that the stability of the commensurate coexistent state is related with the spin 
resonance mode shown in the coexistent phase, as the spin resonance stabilizes the (i-wave 
superconducting state in high- Tc cuprates and CeCoIns. The NMR relaxation rate has also 
been calculated. The calculated result for the system with the coexistent ground state could 
explain well the temperature dependence of NQR data in the coexistent phase of CeRhlns. 
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Appendix A: Mean-Field Equations for the Coexistence of AF and SC 



In the CO phase with A 7^ and m<j 7^ 0, the mean- field Dyson equations in Eq. (fTT]) ^( !T4|l 
involve the four Green's functions, 



G(k,r) 
F"Q(k,r) 

F*Qt(k,r) 



5^^(Tjck+Q.(r)4(0)j) 



(Al) 



The momentm average of F^^(k, r) and F'^^(k, r) = -F'^(k, r)* corresponds to the expectation 
value of the order parameter for the staggered magnetization and the dx^-y^-waye SC gap, 



Norn, = 5^F^Q(k',r = 0), 

k' 

Ak = 5^yk,k'F^(k',r = 0). 



(A2) 



The last term F**^"''(k, r) makes a coupling between F"^(k, r) and F^(k,T) as shown in 
Eq. (m. 

After the Fourier transformation to the imaginary frequency, the generalized Dyson- 
Gorkov equation is written as 



(A3) 



G(k, iun) 




■GW(k,ia;„)" 


F^Q(k,ia;„) 


= 
















_F*Qt(k,ia;„)_ 








with 



A 



1 G^^\k)Ums 
G^^\k + Q)Um, 1 



-G'(o)(-A;)A* 







-GW(-A;-g)A 



G^'\k)A^ 

1 

k+Q -G^'\-k - Q)Um, 



GW(A; + g)Ak+ 
-G^^\-k)Um^ 
1. 



'A 4) 
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where G^'^\k) ^ = iun — and G^'^\k + Q) ^ = iun — £k+Q- The Green's function is given 
by using the cofactor expansion, 

An 



detA 

Al2 



detA 
A 



detA 
Au 
detA 



G(°)(k,i 



lUJr, 



(A5) 



where detA is obtained by using the cofactor Ay (j = 1,4). If the superconducting gap 
function satisfies Ak_|_Q = — hke the case of = A(cos A;^. — cos ky), the secular equation 



\G^'>\k)G^''\k + Q)\-^ detA 



- El 



(A6) 



gives rather simple two eigenvalues E\f_± as 



EL 



4 



+ A^ 



Then, the Green's functions are given as 



(AT) 



k+Q 



] - {iun - ei,+Q){Umsy + (iw„ + £:k)|Ai, 



D(k,iw„)G(k, iun) 

D(k,icu„)F''*^(k,icu„) = (-f/m^)(-ia;„ - e]^){-iujn - ^k+q) + (UrrisY - {-Ums)\Ai, 

D(k, iLOn)F'{k,iUn) 

Z}(k,ia;„)F*Qt(k,ia;„) 



4+Q + (f/"^.)' + |AkP] 



(A8) 



with D(k, iun) - 
the following mean-field equations are obtained 



i^n)' - EL 



. Using the Green's functions in Eq. (1A8I) . 



He 1 \ ^ 



1 



Ek+ , , -E^k+ 
■ tanh ■ 



A, 



E^- . 1 -^k- 

tanh ■ 



2r 2Ek 

-E^k+ , , -E'k+ 
tanh ■ 



2T 
^k 



2E. 



'k,k' 



Ak' , , -Ek' 
— = — tanh 

2£'k'4- 



2T 2Ek'- 



2r 2E, 

Ak' ^ , -Ek'^ 
tanh ■ 



tanh ■ 



E. 



2T 



2T 



where is the momentum summation within the magnetic Brillouin zone. 



(A9) 
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Appendix B: Calculation of spin fluctuaion spectrum in the presence of AF and SC 
long-range order 



In the CO phase, there are two sites in the magnetic unit cell due to the staggered moment 
rris. Then, the area of magnetic Brillouin zone is a half of the original one, and the ordering 
wave vector Q becomes a reciprocal lattice vector. Therefore, momenta of initial and final 
states of relevant spin fluctuation can differ by the reciprocal lattice vector as the Umklapp 
process. Since this fact is independent of components of spin operators, we introduce a 
generalized form of the fluctuation with a, /3 = 0, x, y, z, components, as follows, 

X"^(q,q',ifi„)= / rfre^^"^(T. 6S^{t)6S^_^,{0) ), (Bl) 

Jo L -I 

where = corresponds to the charge operator. Furthermore, when the direction of 
staggered moment is parallel to ^-direction, the fluctuation matrix is decomposed into two 
4x4 matrices formed by (0, z)-components and (x, ?/)-components. 

The transversal components accompanying the contribution of spin wave mode are given 

by 

X^^(q,q',ifi.) = J]x'^-'(q,q',ifi.), (B2) 
X^nq,q',i^n) = 5^aa'r'^'(q,q',in„), (B3) 
X^nq,q',ifin) = i5^a'r'^'(q,q',il^„,), (B4) 
X^"(q,q',ifin) = -iJ2^r^'i^^^'^^n), (B5) 



with 



'(q, q', ifi„) = / dT e'^"-- 5^(T,[cL(r)ck+q_.(r)4,+q,„,,c,q.,]). (B6) 

k,k' 

Through the analytical continuation to the real axis, the transversal susceptibility matrix 
obtained within RPA is given by 

xipM. u + iS) = [i- 2^^(q, cu + iS)U^]-'i^ici, cu + iS), (B7) 
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with the interaction matrix 



f/ 

t/ 

f/ 

t/ 



fB8) 



and the irreducible transversal susceptibility 



X^(q,w + i(5) 



rnq,q + Q) 

rnq + Q,q + Q) rnq + Q,q) o 



x^nq + Q,q + Q) 



q) 





X^"(q,q + Q) 

X^nq+Q,q) 



X^^(q, q) 



where the frequency dependence of every matrix element is omitted in the right hand side. 
The matrix elements of irreducible transversal susceptibility are given with a; = cu + 15, as 
follows, 



x'^'^(q,q,w) = x™(q,q,'^) 



a 



k,k+q{'^k,k+q 



2No ---- ^ ---- E^^ - EkH-q+ + ^ 



+ 



1 1 - /(^k+q+) - /(^k+) ^ 1 - /(^k+q+) - /(^k+ 



k,k+q , 

-C/k+ + -C'k+q^ 
/(Ek+q-) - /(^k-) 



^k 



k+q+ ^ ^ 



^1 C,//- ( 1 - /(^k+q-) - /(^k-) ^ 1 - /(^k+q-) - /(^k^ 



2 '^'^+''\ ^k-+^k+q-+W E^^+E- 



k+q- 



rqI+~ ( /(-^k+q-) - f{El,+ ) , /(-£^k+ q~) - /(-£^k+ 
"'"f-^k,k+qt'-'k,k+q ' ~ ~ + _ 



V ^k 



^1 



k+q- 



E]i+ — -Ek+q- — W 



_ / 1 - /(^k+q-) - /(^k+) ^ 1 - /(^k+q-) - /(^k4 
"^■^k.k+q I ^ ^ ~ ^ 



Ek+ + -E'k+q- + ^ 



-E'k+ + -E'k+q- 



}], (BIO) 
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rnq,q±Q,cD) = -x^^(q,q±Q,(I;) 

^V^r,^- r^/++ /(^k+q+) - /(£'k+) 

^^^0 ^ - ^k+q+ + 

l^n++ ( 1 - /(^k+q+) - /(^k+) 1 - /(^k+q+) - /(^k+) 



2 k,k+q I _^ -Ek+q+ + 1^ -E'k+ + -£'k+q+ 



rpl /(-Ek+q-) - /(-S'k- 

iy//- f 1 - /(^k+q-) - /(^k-) 1 - /(^k+q-) - /(^k^ 



2 k,k+q ^ ^ -Ek+q- + W -Ek- + -Ek+q- ~ ^ 

;^/+_ / /(^k+q-) - /(^k+) /(^k+q^) - /(^k+) 
+7k,k+qt^k,k+q I p p ~ P P 

, rrII+~ f 1 ~ /(-S'k+q-) - f{E]^+) 1 - /(-Ek+q^) - f{El^+) 1-11 X 

+-'k,k+q p , p p MP r, ^ ^ 

with coherence factors due to antiferromagnetism 

± _MlJ- ^Ik^lk+q ± (f^"^) 

•^k.k+q ~ 2 



/3^k^„ = ^ I 1 



1 I , ^Ik^lk+q ± {U' 



m] 



^k,k+q 2 



1 / t/m Um 



m)2 



and coherence factors due to superconductivity 

SS,,, = i f 1 + :5!S^|a:^±^i^) . (B13) 

^ \ -C^k5-C^k+qj7 / 



q/ZC'? _ 1 M _ -^kg-^k+q»7 + AkAk+q 
•^k.k+q - 2 p p 

^ \ -C^k^-C^k+qr? 



(B14) 



-'k.k+q - 2 P P ' '''^+'1 ~ 2 P P ■ ^ ' 

These matrix elements of transversal irreducible susceptibilities are calculated through the 
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irreducible susceptibility of x'^'^'(q, q', if2„), as follows, 



T 



^'^'^'(q,q',ifi„) = -5,,,,^^^[Vq'{G'(k,ia;0G'(k + q,ia;i + i(]„) 



k I 



-F^Q(k, ia;/)F^-Q(k + q + Q, ia;, + il]„)} 
-5q-q'+Q^{G'(k, ia;Oi^^~*^(k + q + Q, icj; + if]„,) 

-F^Q(k, ia;OG'(k + q, icu^ + ifi^} 
+(5q_q.{F^(k, ia;i)F^^(k + q, iui + il]„) 

-F*Q(k, ia;,)F*Qt(k + q, ic^, + ifi,,)} 
+Vq'+Q^{^*^(k, iw«)^''^(k + q + Q, icJi + iVln) 

-F^(k, ic^O^'^^^lk + q - Q, ia;; + il^„)}^B16) 

Similarly, the dynamical susceptibility of charge and z-component of spin are estimated 
within RPA as 



x|p^(q,a; + i5) = [i-2:fll(q,a; + i(5)t/ll] ^^"(q, u; + 15), 
with the interaction matrix 



(B17) 



-u 














-u 














u 














u 



(B18) 



and the irreducible spin-susceptibility 



x"(q,a; + i5) 



X°°(q,q) 



X°^(q,q+Q) 

x°°(q + Q,q + Q) x°^(q + Q,q) o 

r°(q,q + Q) r'(q,q) 

r°(q + Q,q) r^(q + Q,q + Q) 



(319) 
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These matrix elements are given as 



X°°(q, q,w) 



-E 



a 



k,k+q 



^k,k+q 



ql++ 
'-^k,k+q 



+ - 



+ 



+ - 



^k,k+q 
'-^k,k+q 

"-"k.k+q 

"-^k.k+q 

^11 

^k,k+q 

qll 

'-'k,k+q 



-Ek+ — -E'k+q+ 
1 - /(i?k+q+) - /(^k+) , 



-E'k+ + -E'k+q+ + ^ 
/(^i^k+q-) - /(^k-) 

Ek- — Ek+ci- + ^ 

1 - /(^k+q-) - /(^k 



k+q 



+ W 



+/^ltk+q{ 



"-"k.k+q 
"-^k.k+q 



/ /(^^k+q-) - /(^k+) 
V -E'k+ — -Ek+q- + 



+ 



'-'k.k+q 

qll+- 
"-^k.k+q 



1 - /(^k+q+) - /(^k+) 



k+ 



k+q+ 



1 - /(^k+q^) - /(^k^ 
-E'k- + -£^k+q- — ^ 

, /(^k+q-) - /(^k+ 



k+ 



E, 



k+q- 



/(Ek+q_) - ^ 1 - /(^k+q^ 



k+J 



-E'k+ + -E] 



k+q- 



UJ 



E, 



k+ 



^1 



k+q- — 1^ 



}], (B20) 
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X°^(q,q±Q,cD) 
X^°(q,q± Q,^) 



2N, 



1 

+ 2 



k+q+. 



■k+j 



^k,k+q 

C//++ 
"-^k.k+q 

'-'k.k+q 

'-'k.k+q 

^11 

•-^k.k+q 

qll 

"-"k.k+q 



-E'k+ — -E'k+q+ + ^ 
1 - /(Ek+q+) - /(^k+) , l-/(^k+q+) 



k+J 



7k,k+q{ 



-E'kH- + -E'k+q+ + 1^ 

/(i^k+q-) - /(^k-) 
-Ek- — -Ek+q- + 

1 - /(^k+q-) - /(^k- 

-Ek- + -E'k+q- + a) 

/(Ek+q_) - /(^k 



^k 



k+q+ 



^k,k+q 
•-^k.k+q 



+ 



rii+- 

^k,k+q 
"^k.k+q 



-E'k+ — -E'k+q- - 

/(i^^k+q-) - /(^k^ 




-E'k+ — -Ek+q- — 
/(^k+q-) - /(^k+J 



^1 



k+ 



^1 



k+q- 



E, 



k+q- 



UJ 



and coherence factors Cl.^^^^ and C(/^^^ 



•-^k.k+q 



^k,k+q 




-E'k^-E'k+qr; " AfcAk+q 

-E'k^-E'k+q?? 
-E'kg-E'k+qr? ~ AkAk+q 



}], (B21) 



(B22) 



(B23) 



-E'k^-E'k+qry 

In the absence of the magnetic long-range order (m = 0), the coherent factor 7kk+q 
in Eq. (IB12P vanishes, which makes the off-diagonal components of the susceptibility 
X°^(q, q ± Q, if2„) with a ^ P zero and only the diagonal component x""" {q^, q^, iQn) with 
a = X, ?/, 2,0 remain finite. In addition, another coherent factors /^kk+q ^^^'^ vanishes to 
drop the last four terms from the diagonal components x"°(ct, q, if^n) with a = x,y, z,0 in 



Eq. flBlOp and Eq. flB20p . Further it makes the other coherent factor Q^kk+q 



a 



k,k+q 



which recovers the rotational symmetry in spin space and our results of Eq. flBlOp and 
Eq. (]B20I) reduces to the symmetric form of the pure superconducting case^. 
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